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We present a search for Higgs boson in final states with two oppositely charged leptons and large 
missing transverse energy as expected in H — > WW — > Ivl'v' decays. The events are selected 
from the full Run II data sample of 9.7 fb _1 of pp collisions collected with the DO detector at the 
Fermilab Tevatron Collider at y/s = 1.96 TeV. To validate our search methodology, we measure the 
non-resonant WW production cross section and find aww = H-6 ± 0.7 pb, in agreement with the 
standard model prediction. In the Higgs boson search, no significant excess above the background 
expectation is observed. Upper limits at the 95% confidence level on the Higgs boson production 
cross section are therefore derived. Within the standard model, the Higgs boson mass range 159 < 
Mr < 176 GeV is excluded while the expected exclusion sensitivity is 156 < Mh < 172 GeV. 
For a mass hypothesis of Mh = 125 GeV, we exclude Higgs boson production cross sections 4.1 
times larger than the standard model expectation, which is compatible with the presence of a 
Higgs boson at this mass. Within a theoretical framework with a fourth generation of fermions, 
the mass range 125 < Mh < 218 GeV is excluded. The search results are also interpreted in the 
context of fermiophobic Higgs boson couplings, which yields an exclusion of fermiophobic Higgs 
boson production cross sections 3.1 times larger than the expectation for Mh = 125 GeV. 

PACS numbers: 14.80. Bn, 13.85.Qk, 13.85.Rm, 14.65. Jk, 14.80.Ec, 14.70.Fm 



INTRODUCTION 

Spontaneous breaking of the SU{2) x U(l) electroweak 
symmetry explains why the W and Z weak vector bosons 
are massive particles. However the details of the sym- 
metry breaking mechanism are yet to be fully explored. 
In the standard model (SM), it results from the exis- 
tence of a single elementary scalar field doublet that ac- 
quires a non-zero vacuum expectation value. After ac- 
counting for the mass of the weak vector bosons, one 
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degree of freedom remains, manifesting itself as a single 
scalar particle, the Higgs boson. Its mass, Mh, is a free 
parameter of the model. A lower limit of 114.4 GeV 
was set on Mh by the CERN LEP experiments [lj]. 
This experimental constraint was extended by the com- 
bined results from the CDF and DO experiments that 
excluded the Higgs boson mass range from 156 GeV to 
177 GeV Upper (lower) limits of 131 (122) GeV [| 

and 128 (121.5) GeV [5J have then been established by 
the ATLAS and CMS Collaborations, respectively. These 
exclusion limits and those reported hereafter are all de- 
fined at the 95% C.L. In both Ref. H and Ref. @, ex- 
cesses above background expectations at the five stan- 
dard deviation (s.d.) level have been reported, consistent 
with the observation of a Higgs boson of Mh ~ 125 GeV. 
The CDF and DO Collaborations have reported excesses 
above background expectations in the H — > bb search 
channels [H, Q- Their combination yields an excess at 
the three s.d. level, consistent with the production of a 
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Higgs boson of mass M H « 125 GeV [8(. 

In this Letter, we present a search for the SM Higgs bo- 
son in final states containing two oppositely charged lep- 
tons {££' = efj,, ee, or /z//, where small contributions from 
lcptonic t decays are also included) and missing trans- 
verse energy ($r). The search relies on the full Run II 
data set of 9.7 fb _1 of pp collisions collected with the DO 
detector at the Fermilab Tevatron Collider at y/s = 1.96 
TeV. This analysis supersedes our previously published 
results in the same final states, obtained after analyzing 
5.4 fb" 1 @ and 8.6 fb" 1 of integrated luminosity. 
A similar search has been conducted by the CDF Col- 
laboration using 4.8 fb -1 of integrated luminosity fill ]. 
The results from Refs. 0, [ll| have been combined in 
Ref. • More recently, searches in dilepton plus missing 
transverse energy final states have been conducted by the 
ATLAS [3 and CMS [H Collaborations using 4.7 ftr 1 
and 4.6 fb -1 of integrated luminosity, respectively. 

The main Higgs boson production and decay channel 
resulting in opposite charge dilepton plus Iftr final states 
at the Tevatron is the gluon fusion production, gg — > H, 
with subsequent decay H — > W + W~ — > £ + v£'~v' , where 
one of the W bosons is virtual for Mh < 160 GeV. 
This final state receives additional contributions from 
Higgs boson production via vector boson fusion (VBF), 
qq' — > qq'VV — > qq' ' H, and from production in associa- 
tion with a vector boson, qq' — > VH (V = W,Z). The 
primary source of background is due to diboson produc- 
tion, in particular the non-resonant pp — > WW process. 
Other background sources are the Drell-Yan (DY) pro- 
cess, pp — > Z/^i* — > 11, with a mismeasured IpT, the 
lcptonic decays of top-quark pairs (£i), M^+jets/7 and 
multijet events in which jets (photons) are misidentificd 
as leptons (electrons). 

The initial selection of Higgs boson candidate events 
is based on the reconstruction of two high transverse 
momentum (pt) leptons. This selection is followed by 
additional requirements, involving and the usage of 
multivariate techniques based on boosted decision trees 
(BDT) (3], to suppress the large DY background. To 
increase the sensitivity, the events are separately ana- 
lyzed according to the lepton flavors (ee, efi, and h/j,) 
and jet multiplicity, and they are also categorized into 
WW-enriched and VyW-depleted sub-samples. Addi- 
tional BDTs are trained to separate the signal from the 
remaining background events. To demonstrate the va- 
lidity of the techniques used in this search, we use sim- 
ilar BDTs to measure the cross section for the SM non- 
resonant WW production cross section. For the Higgs 
boson searches, the outputs of the BDTs are the final 
discriminants used for the statistical interpretation of the 
data, within the SM framework, but also in the contexts 
of a fourth generation of fermions and a fermiophobic 
Higgs boson. These models are described in more detail 
in later sections. 



DETECTOR AND OBJECT RECONSTRUCTION 



The DO detector used for Run II (2002 - 2011) is de- 



scribed in detail in Ref. [151 ] - The innermost part of the 



detector is composed of a central tracking system with 
a silicon microstrip tracker (SMT) and a central fiber 
tracker embedded within a 2 T solenoidal magnet. The 
tracking system is surrounded by a central preshower de- 
tector and a liquid- argon/uranium calorimeter with elec- 
tromagnetic, fine, and coarse hadronic sections. The 
central calorimeter (CC) covers pseudorapidity [HI \rj\ 
< 1.1. Two end calorimeters (EC) extend the coverage 
to 1.4 < |?y| < 4.2. The pseudorapidity gap between the 
ECs and CC is covered by scintillating tiles. A muon 
spectrometer, with pseudorapidity coverage of |?y| < 2, 
resides outside the calorimetry and is comprised of drift 
tubes, scintillation counters, and toroidal magnets. Trig- 
ger decisions are based on information from the tracking 
detectors, calorimeters, and muon spectrometer. 

Electrons are reconstructed as isolated clusters in the 
electromagnetic calorimeter, and required to spatially 
match a track in the central tracking system. They have 
to pass a BDT (ee channel) or likelihood (e/i channel) 
criterion (collectively called electron quality later in the 
text) that accounts for calorimeter shower shape observ- 
ables, a spatial track match probability estimate, and the 
ratio of the electron cluster energy to track momentum 
{E/p). Electrons are required to be in the acceptance of 
the calorimeter (|?y| < 1.1 or 1.5 < |?7| < 2.5). In the di- 
electron channel, events with one electron in the EC are 
treated separately from events with both electrons in the 
CC. Events where both electrons arc in the EC arc not 
considered due to a large background and a small signal 
contribution. 

Muons are identified by the presence of at least one 
track segment reconstructed in the acceptance (I77I < 2.0) 
of the muon spectrometer, that is spatially consistent 
with a track in the central tracking detector. The mo- 
mentum and charge are measured by the curvature in 
the central tracking system. To select isolated muons, 
criteria based on the momenta of central tracks emitted 
in the approximately same direction as the muon and 
criteria based on the energy deposited around the muon 
trajectory in the calorimeter are employed. The number 
of hits in the wire chambers and in the scintillators are 
combined to define a muon quality variable used in the 
final stage of the analysis. 

Jets are reconstructed from energy deposits in the 
calorimeter using an iterative midpoint cone algo- 
rithm 17 1 with a cone radius 1Z = 0.5 [Hfl . The jet 
energies are calibrated using transverse momentum bal- 
ance in 7+jet events [3]. Jets are considered in this 
analysis only if they have pr > 20 GeVand |?7| < 2.4. 
Each jet is also required to be matched to at least two 
tracks associated to the pp interaction vertex. 
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The IpT and its direction are obtained from the vector 
sum of the transverse components of energy deposits in 
the calorimeter, corrected for the differences in detector 
response of the reconstructed muons, electrons, and jets. 



DATA AND SIMULATED SAMPLES 

Signal and SM background processes except multijet 
are simulated with pythia [20] or ALPGEN [21j genera- 
tors, with PYTHIA providing showering and hadronization 
in the latter case, using the CTEQ6L1 [22( parton distri- 
bution functions (PDFs) , followed by a detailed GEANT3- 
based [23| simulation of the DO detector. In order to 
model the effects of multiple pp interactions, the Monte- 
Carlo (MC) samples are overlaid with events from ran- 
dom pp collisions with the same luminosity distribution 
as data. Then, these events are reconstructed with the 
same software as used for the data. Jet energy calibra- 
tion and calorimeter response to unclustered objects are 
adjusted in simulated events to match those measured in 
data. Corrections for residual differences between data 
and simulation are applied to electrons, muons, and jets 
for both identification efficiencies and energy resolutions. 

Higgs boson signal samples are simulated using PYTHIA 
for 100 < M H < 200 GcV in increments of 5 GeV, and 
for 200 < M H < 300 GeV in increments of 10 GeV. 
For gg — ► H production, the cross section is calculated 
at next-to-next-to-leading order with resummed next-to- 
next-to-leading logarithm (NNLO+NNLL) [H], for VBF 
at NNLO [iH, and for VH at NNLO [H]. All signal cross 
sections are^ computed using the MSTW 2008 NNLO 
PDF set " 



27] . The distribution of the Higgs boson px 



for gg H process is weighted to match the calculation 
of the hqt generator, which has NNLO and NNLL accu- 
racy (28|. The Higgs boson branching ratio predictions 
are taken from HDECAY [2§| . 

The IV+jets and Z+jets backgrounds are modeled us- 
ing alpgen. The W+jets and Z+jets processes are nor- 
malized using the NNLO cross section calculations of 
Ref. O ■ The px distribution of Z bosons is weighted to 

taking into 



match the distribution observed in data 31 



account its dependence on the number of reconstructed 
jets. The px distribution of W bosons is weighted to 
match the measured Z boson px spectrum, corrected for 
the differences between the W and Z px spectra pre- 
dicted in NNLO QCD |32|. In the ee and e/z chan- 
nels, the FV+jets simulation includes contributions from 
events where a misidentified electron originates from a 
jet or a photon. The size of each of these contributions 
is corrected such that the distribution of the number of 
hits in the innermost silicon layer, associated to the elec- 
tron track, matches that observed in a VF+jets enriched 
control sample. 

The it process is modeled using alpgen with shower- 
ing and hadronization provided by pythia, and normal- 



ized to the approximate NNLO cross section [331 ] . 

Diboson production processes {WW , WZ, and ZZ) 
are simulated using pythia, normalized to NLO cross 
sections [341 ]. An additional correction, determined us- 
ing the POWHEG generator [2E], accounts for Z/j* in- 
terference in WZ production. For the irreducible back- 
ground arising from WW production, the px of the dibo- 
son system is modeled using the MC@NLO simulation [3r| . 
and the distribution of the opening angle of the two lep- 
tons is corrected for the contribution of the non-resonant 
gg -> WW process [37 1. 

The background due to multijet production, where jets 
are misidentified as leptons, is determined from data by 
inverting some of the lepton selection criteria. All other 
event selection criteria are applied in order to model the 
kinematic distributions of the multijet background in the 
signal region. In the \x\i channel, the opposite-charge re- 
quirement for muons is reversed, and a correction for the 
presence of non-multijet events in the like-sign sample, 
estimated from simulation, is applied. For the ee and e/i 
channels, the electron quality requirement is reversed, 
and the normalization is determined from control sam- 
ples in which the leptons have the same charge. 



EVENT PRESELECTION 

A first selection is applied on the data by requiring 
two high-pr leptons, that have opposite charge and that 
originate from the same location, within 2 cm, along 
the beamline. In the ee and nfi channels, the highest- 
px and next highest px leptons are required to satisfy 
Px > 15 GeV and p^ > 10 GeV, whereas in the e[i 
channel, p T > 15 GeV and p T > 10 GeV are required. 
Additionally, in the ee and final states, the dilepton 
invariant mass Mu is required to be greater than 15 GeV. 
A veto against additional leptons is applied to prevent 
overlap with dedicated Higgs searches in trilepton final 
states [38|] . These criteria define the "preselection" stage 
of the analysis, and they select samples dominated by 
DY production. Most events selected at this level pass 
single-lepton trigger conditions. But, as no specific trig- 
ger requirement is made, the trigger acceptance with re- 
spect to off-line selections is enhanced to rj 92% for the 
fifj, channel and « 100% for the ee and e/x channels, due 
to additional events passing lepton+jets or dilepton trig- 
gers. The remaining trigger inefficiency is modeled in the 
simulation by corrections derived from Z — > ££ samples 
selected with different trigger requirements. The prese- 
lected samples are further subdivided according to the 
number of jets in the event. Namely 0-, 1-, and (> 2)-jet 
multiplicity bins are considered. Dividing the analysis 
into different jet multiplicity bins significantly increases 
the sensitivity of this search as the signal and background 
compositions are different in each sample. 

To correct for any possible mismodcling of the lepton 
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TABLE I: Observed and expected number of events after preselection in the ee, efj,, and jj.fi final states. The signal is for a Higgs 
boson mass of 125 GeV. The uncertainty quoted on the background combines both statistical and systematic uncertainties, 
after the normalization procedure described in the main text. 



Data Total background Signal Z/j* —> ee Z/j* —> /i/i Z/j* —> tt tt VU+jets Diboson Multijet 



ee: 


659570 


664460 


± 


13290 


16.1 


653263 




5494 


210 


795 


945 


3752 


e\i\ 


14936 


15142 


± 


303 


16.6 


408 


1211 


8671 


537 


1225 


906 


2184 


Mi: 


811549 


818269 


± 


16370 


18.7 




807642 


6459 


356 


438 


1314 


2060 



reconstruction and trigger efficiencies, and to reduce the 
impact of the luminosity uncertainty, scale factors are 
applied to the MC samples at the preselection stage to 
match the data. The Z boson mass peak regions in the 
preselected samples are used to determine normalization 
factors. Their differences from unity are found to be con- 
sistent with the luminosity uncertainty of 6.1 % [5^] . This 
procedure is repeated for each jet multiplicity to derive 
jet-bin-dependent DY background normalizations to cor- 
rect for possible mismodeling of the DY jet multiplicity. 

The number of events after the preselection is pre- 
sented in Table U [!(| ; all sub-samples are dominated by 
DY production. Figure DTa) shows the dilepton invari- 
ant mass distributions [40j for data and the background 
prediction for the combined sub-samples. 



ANALYSIS USING DECISION TREES 

In the ee and fifi channels, BDTs are trained for each 
Higgs boson mass value and each jet multiplicity bin 
to discriminate between the signal and the dominant 
DY background. The input variables to these "anti-DY 
BDTs" arc kinematic quantities, such as the lepton mo- 
menta, the azimuthal opening angle between the two lep- 
tons, IpT, variables that take into account both Tpr and 
its direction relative to a lepton or a jet, and observables 
that differentiate between real and misreconstructcd JpT- 
This multivariate technique follows the method defined 
in the previous publication [l(J where more details on the 
BDTs' input variables are given. The final selection stage 
for the ee and fifi channels is obtained by applying cuts 
on the anti-DY BDT discriminants [4(J. The thresholds 
are chosen to obtain similar background rejection as the 
cut-based rejection employed in Ref. 

In the e\x channel the final selection stage requires 
Mf in > 20 GeV and M T2 > 15 GeV [4(J, where 
M™ m is the minimum value, over the two possible 
lepton choices, of the transverse mass, Mt(£,$t) = 

^2 ■ ■ $t ■ [1 — cos A<p(£, ]fcr)] , an d Mt2 is an exten- 
sion of the transverse mass suitable for final states with 



two visible and two invisible particles [41 



The number of events at the final selection stage can 
be found in Table UH and the distribution of the angular 
separation between the leptons, combined for all dilepton 



final states, TZ(£ + £~), can be seen in Fig.^b) 0. 

At the final selection stage, a scries of new BDTs is 
built: the "WW-BDTs" are trained to separate the non- 
resonant WW production from the other backgrounds, 
while the "final BDTs" are trained to separate the signal 
from all the backgrounds. In the former case, the Higgs 
signal is not used in the training. These BDTs rely on the 
same input variables as for the anti-DY BDTs, but sup- 
plementary variables are added characterizing the lepton 
reconstruction quality and the lepton isolation, to dis- 
criminate against the instrumental backgrounds (multi- 
jet and W+jets backgrounds). Outputs of jet 6-tagging 
multivariate discriminants [42| are also added as inputs 
to separate the signal from the tt background. 

Using the JUVU-BDT discriminants, we split the 0- and 
1-jet samples into VUVU-depleted and TUVU-enriched re- 
gions for the ee and fifi analyses. In the e\i channel, 
splitting only the 0-jet sample according to the lepton 
reconstruction quality achieves a sufficiently pure sepa- 
ration of the data sample into a lUlU-depleted and WW- 
enriched sub-samples. The final BDTs are then trained 
separately for each jet multiplicity bin, for each dilepton 
final state, and for the VUVU-dcplcted and VUVU-cnrichcd 
samples, resulting in 14 BDTs for each mass hypothe- 
sis [40(. The outputs of these BDTs are used as final 
discriminants. Figure [TJc) shows the BDT distributions 
of the 14 sub-samples summed in bins with similar signal 
to background ratios (s/b). 



SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are estimated for each final 
state, background, and signal process. They can affect 
only the normalization or both the normalization and the 
shape of the final discriminants. 

Sources of systematic uncertainty that affect only the 
normalization arise from the overall normalization un- 
certainty due to theoretical inclusive cross sections of 
Z+jcts (4%), W+jcts (6%), diboson (6%) and tt (7%) 
processes; multijet normalization (30%); the VF+jets jet- 
bin-dependent normalization (15%— 30%); the Z+jets jet- 
bin-dependent normalization (2%-15%); and the model- 
ing of the IfjT measurement for the Z+jets background 
(5%-19%). 

The uncertainties on a(gg — > H) production are esti- 
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FIG. 1: [color online] The (a) dilepton invariant mass at preselection level, (b) the angular separation 1Z(£,£) between the 
leptons at preselection level, and (c) output of the final BDT discriminants after the final selection, sorted as a function of 
signal over background ratio. In (a), the last bin includes all events above the upper bound of the histogram. In these plots, 
the hatched bands show the total systematic uncertainty on the background predictions, and the signal distributions are those 
expected from a Higgs boson of mass Mh = 125 GeV. It is scaled by a factor 20 in (c). 



TABLE II: Expected and observed number of events after the final selection in the ee, e/x, and fifj, final states. The signal is for 
a Higgs boson mass of 125 GeV. The numbers in parentheses correspond to the efficiency of the final selection with respect to 
the preselection, shown in Table U for both the total background and signal. The uncertainty quoted on the total background 
combines both statistical and systematic uncertainties. 
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Total backj 


•round 


Si 


gnal 


Z — > ee 


Z — > /.t/i 


Z — v tt 


tt 


VF+jcts 


Diboson 


Multijet 
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1882 


1859 


± 


205 


(0.3%) 


7.5 


(46.8%) 


746 




55 


151 


518 


371 


18 


Ojet 


1289 


1317 


± 


145 


(0.2%) 


4.6 


(64.8%) 


528 




32 


12 


424 


307 


13 


1 jet 


379 


343 
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38 


(0.4%) 
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152 
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199 
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22 


(1.7%) 
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(27.5%) 


65 




16 
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1954 


1960 
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19 
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mated following the prescription described in Ref. [43(, 
i.e., by considering as uncorrelated the scale uncertain- 
ties of the NNLL inclusive 0, Q, NLO > 1 jet 0, 
and NLO > 2 jets 46j cross sections. This prescription 
results in the following covariance matrix for the the ex- 
clusive production of gg — > H + jet, +1 jet, and +2 jets 
or more , respectively: 
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H + jet, +1 jet, and +2 jets or more, respectively. 



(26.6%) 2 
-(28.3%) 2 




-(28.3%) 2 
(41.8%) 2 
-(20.5%) 2 





(20.5%) 2 
(33.0%) 2 



(1) 



ob- 



Thc PDF uncertainties for gg — > H production, 
taincd using the prescription from Refs. 
7.6%, 13.8%, and 29.7% for the exclusive production of 



24, W5J, are 



We also consider sources of systematic uncertainty 
that affect the shape of the final discriminant distribu- 
tion (and we quote here the average fractional uncer- 
tainty across bins of the final discriminant distribution 
for all backgrounds): jet energy scale (4%); jet resolu- 
tion (0.5%); jet identification (2%); jet association to 
the hard-scatter primary pp interaction vertex (2%); b- 
tagging (< 2%); and W+jets modeling (10%-30%), de- 
pending on jet multiplicity bin and final state. The sys- 
tematic uncertainties due to the modeling of pr(WW) 
and Acf> between leptons, and the pr of the vector boson 
from the y+jets production, are at the level of < 1% and 
taken into account. 
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MEASUREMENT OF THE NON-RESONANT 
pp -> WW CROSS SECTION 

To validate the analysis techniques employed to search 
for the Higgs boson, a measurement of the non-resonant 
WW production cross section is performed. This is mo- 
tivated by the fact that WW production yields similar 
particle content and topology as the Higgs boson signal. 
The same analysis methods are employed as for the Higgs 
bosons searches, and same sources of systematic uncer- 
tainty are accounted for, but the outputs of the WW 
discriminants, described in the "Analysis using decision 
trees" section, are considered. The WW cross section is 
obtained as the result of a maximum likelihood fit to the 
data, with maximization over the WW signal normal- 
ization and over the systematic uncertainties treated as 
nuisance parameters, as for the SM Higgs boson search 
results described in the next section. The measurement 
is carried out using discriminants from the three dilepton 
final states, in the 0- and Tjet multiplicity bins. Figure^] 
shows the combined output distribution of these discrim- 
inants rebinned according to s/b and after the ex- 
pected backgrounds have been subtracted. The measured 
value a W w = 11.6 ± 0.4 (stat) ± 0.6 (syst) pb is in 
agreement with the SM prediction of 11.3 ± 0.7 pb [34| . 
The possible presence of a SM Higgs boson of 125 GeV 
in the data is not accounted for, but it is expected to bias 
this measurement upward by ~ 0.1 pb. 



V> I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 

| D0, 9.7 fb' 1 , II + E T 
> 400 - -1,- Data - Bkgd 
- ■ Signal (WW) 
300 ±1 s.d. on Bkgd 




■3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 



log 10 (s/b) 

FIG. 2: [color online] The post-fit background-subtracted 
data distribution for the final discriminant, summed in 
bins with similar signal to background ratios, for the WW 
cross section measurement. The uncertainties shown on the 
background-subtracted data points are the square roots of the 
post-fit predictions for signal plus background events in each 
bin, representing the expected statistical uncertainty on the 
data points. Also shown is the ±1 s.d. band on the total 
background after fitting. 



SM HIGGS BOSON SEARCH RESULTS 

Table iHland Fig. [Tic) demonstrate good agreement be- 
tween the data and the predicted background, in both the 
numbers of selected events and the distributions of final 
discriminants. The modified frequentist CL S method (47| 
is employed to set limits on SM Higgs boson production, 
where the test statistic is a log-likelihood ratio (LLR) for 
the background-only and signal+background hypotheses. 
The LLR is obtained by summing the LLR values of 
the bins of the 14 BDT outputs from the different sub- 
channels. In the LLR calculation the signal and back- 
ground rates are functions of the systematic uncertain- 
ties which are taken into account as nuisance parame- 
ters with Gaussian priors. Their degrading effect is re- 
duced by fitting signal and background contributions to 
the data by maximizing the profile likelihood function for 
the background-only and signal+background hypotheses 
separately, appropriately taking into account all correla- 
tions between the systematic uncertainties fiil ]. 

Figure EJa) shows the LLR values as a function of the 
tested Higgs boson mass hypothesis. The LLR values 
expected in the absence of signal and in the presence of 
a SM Higgs boson of mass 125 GeV are also displayed 
for comparison. Figure E^b) [!o| presents expected and 
observed upper limits for a(pp — > H + X) relative to 
SM predictions. For M H = 165 GeV (125 GeV), the ex- 
pected limit is 0.76 (3.4) times the SM prediction and the 
observed limit reaches 0.74 (4.1) in the same units. A SM 
Higgs boson in the mass range 159 < Mh < 176 GeV is 
excluded at the 95% C.L. while the expected exclusion 
sensitivity is 156 < Mh < 172 GeV. In these figures, 
a slight excess of signal-like candidates yields a limit 
roughly one s.d. above the background expectation, in 
the mass range 100 < Mh < 145 GeV. FigureEJc) shows 
a comparison of the BDT output distributions, sorted as 
a function of signal over background ratio, expected for 
the signal of Mh = 125 GeV, and observed in the data 
after subtracting the fitted backgrounds. 

UPPER LIMIT ON gg^t H -> WW AND FOURTH 
GENERATION FERMION INTERPRETATION 

Additional generations of fcrmions can occur natu- 
rally in models of grand unification, CP violation, gauge- 
mediated supersymmetry breaking, and others. Mea- 
surements of the Z boson decay width exclude mod- 
els in which the fourth neutrino mass eigenstate is lighter 
than 45 GeV, but fourth generation models can still 
be accommodated for a large fourth-generation neutrino 
mass. Production of gg — > H occurs via top-quark loops 
in the SM. With respect to the SM, the quarks from 
the fourth generation will provide additional contribu- 
tions to the quark loop diagram, enhancing production 
by a factor of 7 to 9, depending on their masses and 
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FIG. 3: [color online] (a) The observed LLR as a function of Mh- Also shown are the expected LLRs for the background-only 
hypothesis, for the signal+background hypothesis, and the expectation in the presence of a signal of Mh = 125 GeV. (b) 
Excluded cross section, cr(pp — > H + X), at the 95% C.L. in units of the SM cross section as a function of Mh- In (a) and 
(b), the green and yellow shaded bands indicate ±1 and ±2 s.d. uncertainties of the expected observation for the background- 
only hypothesis, respectively, (c) Background-subtracted data distribution for the final discriminants, summed in bins with 
similar signal to background ratios, for Mh = 125 GeV. The uncertainties shown on the background-subtracted data points are 
the square roots of the post-fit background predictions of number of events in each bin, representing the expected statistical 
uncertainty on the data points. Also shown is the ±1 s.d. band on the total background after fitting. 



the Higgs boson mass [5CH52| . A previous combined 
DO and CDF result using up to 5.4 fb -1 of data ex- 
cluded the existence of a SM-likc Higgs boson in the 
mass range between 131 GeV and 204 GeV [H^, as- 
suming the presence of a fourth sequential generation 
of fermions with large masses. Similar searches have 
been conducted by the ATLAS and CMS Collaborations, 
yielding and exclusion of 140 < run < 185 GeV 54| and 
144 < m,H < 207 GeV [H], respectively. 

To test such models, we derive upper limits on the 
gg — >• H —> WW production cross section. The same 
analysis as described in the previous sections is per- 
formed, but the VBF and VH contributions are excluded 
from the overall signal yield when constructing the LLR. 
The upper limits are reported in Fig. 0] [Io| , compared 
to the expected yield of the gg — > H production in two 
models of fourth generation fermions. In the "low-mass" 
scenario, the masses of the fourth generation charged lep- 
ton and neutrino are assumed to be respectively = 
100 GeV and m U 4 = 80 GeV, just beyond the experi- 
mental limits, which yields a reduction by up to 15% in 
the branching ratio for H —> WW. On the contrary, in 
the "high-mass" scenario, where mn = m v4 = 1 TeV, the 
lcptons are too heavy to contribute to the Higgs boson de- 
cay width and the branching ratio for H — > WW remains 
basically unchanged relative to the SM branching ratio. 
For both scenarios, the masses of the fourth-generation 
down-type (m^) and up-type (m^) quarks are fixed to 
m M = 400 GeV and m u4 = 450 GeV plUH- From this 
figure, we derive exclusion of the Higgs boson mass range 
125 < M H < 218 GeV and 125 < M H < 228 GeV, in the 
low-mass and high-mass scenarios, respectively. 
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FIG. 4: [color online] Excluded cross section a(gg — > H) x 
BR(H —¥ WW) in pb as a function of Mh using all chan- 
nels. The red and blue lines correspond to the theoretical 
prediction for a sequential fourth generation assumption in 
the "low mass" and "high mass" scenarios, respectively (see 
main text). The hatched region corresponds to the PDF and 
scale uncertainties on the fourth generation "low mass" pre- 
diction [jsTJ]. The green and yellow shaded bands indicate ±1 
and ±2 s.d. uncertainties of the background-only hypothesis, 
respectively. 



FERMIOPHOBIC HIGGS BOSON 
INTERPRETATION 



The mechanism of electroweak symmetry breaking 
may offer a richer phenomenology than expected in the 
SM. Several Higgs bosons may exist, or the Higgs bo- 
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son(s) may have couplings different from those predicted 
by the SM. In this section, we explore the possibility that 
the lightest Higgs boson does not couple to fermions at 
the tree level, but still behaves like the SM Higgs boson 
for its other properties, in particular for the coupling to 
vector bosons. In this model, the VBF and VH pro- 
duction have the same cross sections as in the SM. The 
main consequences of the vanishing fcrmion couplings are 
the suppression of production via gluon fusion, and the 
enhancement of the branching ratios to vector bosons, 
H — > WW, H — > ZZ, and H — > 77, particularly sizeable 
below the WW threshold, vtlh < 160 GcV. To provide 
masses to the fermions, additional degrees of freedom 
must exist in the Higgs sector, as predicted in models 
with Higgs doublets or triplets [56}, but it is assumed 
that those other particles do not have phenomenolog- 
ical impact in our search. In this model, the CMS and 
ATLAS collaborations exclude a Higgs boson in the mass 

, while masses be- 



range 110 < m H < 194 GeV [57 



low 110 GeV are excluded by LEP experiments |59h6 



and Tevatron experiments |63j, |64( . The same analysis 



steps are performed as described for the SM Higgs boson 
searches, but the various BDTs are retrained, accounting 
for the fermiophobic Higgs branching ratios, computed 
using hdecay, the VBF and VH production at the SM 
rate, and the suppression of gg — ¥ H production. The 
data are in good agreement with background expectation 
and upper limits on the fermiophobic Higgs are derived, 
following the same method as for the SM Higgs. They 
are reported in Fig. [5] [4o| . We obtain a cross section 
upper limit of 3.1 times the fermiophobic Higgs boson 
production cross section for Mjj — 125 GeV, while the 
expected sensitivity is 2.5. 
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FIG. 5: [color online] Excluded cross section, cr(pp — > H + X), 
as a function of Mh using all channels, in units of the Higgs 
boson production rate expected from the fermiophobic Higgs 
boson model (FHM) described in the main text. 



CONCLUSIONS 

We have performed a search for Higgs boson produc- 
tion using final states with two oppositely charged lep- 
tons and large missing transverse energy in the e^, ee, 
and fifj. channels. To validate our search methodology we 
have measured the non-resonant WW production cross 
section, which yields aww = H-6 ± 0.7 pb, in good 
agreement with the SM prediction of 11.3 ± 0.7 pb. For 
the Higgs boson searches, we observe agreement between 
data and the expected backgrounds. We set upper lim- 
its on SM Higgs boson production at the 95% C.L. that 
exclude the mass range 159 < Mh < 176 GeV, while the 
expected exclusion sensitivity is 156 < Mh < 172 GeV. 
For a mass hypothesis of Mh = 125 GeV, we exclude 
4.1 times the expected SM Higgs boson production cross 
section, while the expected sensitivity is 3.4. This upper 
limit is compatible with the presence of a SM Higgs bo- 
son of Mh = 125 GeV. We also interpret our search 
results as cross section upper limit for gg — > H pro- 
duction, which allows us to exclude the mass range 
125 < Mh < 218 GeV in the context of a fourth genera- 
tion of fermions. The search results are also interpreted 
in the framework of a fermiophobic Higgs boson, which 
yields an exclusion of 3.1 times the fermiophobic Higgs 
boson production rate for Mh = 125 GeV, while the 
expected sensitivity is 2.5. 
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TABLE III: Expected and observed numbers of events after the preselection in the ee, efi, and fifi final states in the different 
jet multiplicity bins. The signal is for a Higgs boson mass of 125 GeV. 





Data 


Total background 
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Z -> ee 


Z — > fifi 
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W+jets 


Diboson 


Multijct 
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Ojet 
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51 
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180 
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Ojet 
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> 2 jets 
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78 
16 
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313 
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Event yields after the preselection are shown in Table ITTll The distributions of the dilepton invariant mass (Fig. [6]), 
angular separation between the leptons (Fig. [7]) and missing transverse energy (Fig. [SJ show the good agreement 
between the data and the simulation after the preselection in the ee, efi and fifi channels. Figures [91 and HU1 show the 
distributions of the Drell-Yan (DY) BDT discriminants for the ee and fifi channels respectively. Figure [TT] shows the 
distributions of the variables M™ ln and Mt2, used by the efi channel to reject the DY background. The distributions of 
the dilepton invariant mass (Fig. I12[) . angular separation between the leptons (Fig. I13[) and missing transverse energy 
(Fig. I14[) arc shown for the ee, efi and fifi channels after the final selection. The final discriminant distributions for 
the three channels are shown in Figs. [T~5l [TBI andfTTl Figure [T51 shows the background-substracted data distributions of 
the final discriminants for Mji = 125 GeV and Mh = 165 GeV. Tables HVl fVl and I VII give the expected and observed 
upper limits at the 95% C.L. for Higgs boson production in the SM, in models with a fourth generation of fcrmions 
and assuming fermiophobic couplings, respectively. Figures [T9j [20] and [2TJ show the WW discriminant distributions 
for the ee, efi and fifi channels. Figure [22] shows the distribution of the combination of WW discriminants for the ee, 
e/i and fifi channels. The results of the pp — > WW cross section measurement for the ee, ey y,[i channels, and their 
combination are shown in Fig. 1231 
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FIG. 6: Distributions of the dilepton invariant mass for the (a) ee channel, (b) e/i channel, and (c) \x\i channel after the 
preselection. In (a), (b) and (c) the last bin includes all events above the upper bound of the histogram. In these plots, the 
hatched bands show the total systematic uncertainty on the background predictions, and the signal distributions are those 
expected from a Higgs boson of mass Mh = 125 GeV. 
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FIG. 7: Distributions of the angular separation 1Z(£,£) between the leptons for the (a) ee channel, (b) e/i channel, and (c) 
/i/i channel after the preselection. In these plots, the hatched bands show the total systematic uncertainty on the background 
predictions, and the signal distributions are those expected from a Higgs boson of mass Mh = 125 GeV. 
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FIG. 8: Distributions the missing transverse energy for the (a) ee channel, (b) efi channel, and (c) \i\x channel after the 
preselection. In (a) , (b) and (c) the last bin includes all events above the upper bound of the histogram. In these plots, the 
hatched bands show the total systematic uncertainty on the background predictions, and the signal distributions are those 
expected from a Higgs boson of mass Mh = 125 GeV. 
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FIG. 9: Distributions of the Drell-Yan (DY) BDT discriminant for ee channel in the (a) 0-jet bin, (b) 1-jet bin, and (c) > 
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FIG. 10: Distributions of the DY BDT discriminant for channel in the (a) 0-jet bin, (b) 1-jet bin, and (c) > 2-jets bin. 
The BDTs are trained for a Higgs mass of 125 GeV. The signal distributions are those expected from a Higgs boson of mass 
Mh = 125 GeV. 
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FIG. 11: Af™ ln distribution for the e\i channel in the (a) 0-jet bin, (b) 1-jet bin, and (c) > 2-jets bin. Mti distribution for 
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FIG. 12: Distributions of the dilepton invariant mass for the (a) ee channel, (b) e/x channel, and (c) /xp channel after the 
final selection. In (a), (b) and (c) the last bin includes all events above the upper bound of the histogram. In these plots, 
the hatched bands show the total systematic uncertainty on the background predictions, and the signal distributions are those 
expected from a Higgs boson of mass Mh = 125 GeV. 
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FIG. 13: Distributions of the angular separation 1Z(£ , t) between the leptons for the (a) ee channel, (b) e\x channel, and (c) 
channel after the final selection. In these plots, the hatched bands show the total systematic uncertainty on the background 
predictions, and the signal distributions are those expected from a Higgs boson of mass Mh — 125 GeV. 
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FIG. 14: Distributions of the missing transverse energy for the (a) ee channel, (b) e\x channel, and (c) \i\x channel after the 
final selection. In (a), (b) and (c) the last bin includes all events above the upper bound of the histogram. In these plots, 
the hatched bands show the total systematic uncertainty on the background predictions, and the signal distributions are those 
expected from a Higgs boson of mass Mh = 125 GeV. 
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FIG. 15: Distributions of the final BDT discriminant for Mh = 125 GeV for the ee channel with (a) 0-jet VKVK-depleted, (b) 
0-jet WW-enriched, (c) 1-jet WW-depleted, (d) 1-jet W W-enriched, and (e) > 2-jets. 
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FIG. 16: Distributions of the final BDT discriminant for Mh = 125 GeV for the e/x channel with (a) 0-jet WW-depleted, (b) 
0-jet VKW-enriched, (c) 1-jet , and (d) > 2-jets. 
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FIG. 17: Distributions of the final BDT discriminant for Mh = 125 GeV for the channel with (a) 0-jet WW-depleted, (b) 
0-jet WW-enriched, (c) 1-jet WW-depleted, (d) 1-jet WW-enriched, and (e) > 2-jets. 




FIG. 18: The background-subtracted data distributions for the final discriminants, summed in bins with similar signal to 
background ratios, for (a) Mh = 125 GeV, and (b) Mh = 165 GeV. The uncertainties shown on the background-subtracted 
data points are the square roots of the post-fit background number of events predictions in each bin, representing the expected 
statistical uncertainty on the data. Also shown is the ±1 s.d. band on the total background after fitting. 
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TABLE IV: Expected and observed upper limits at the 95% C.L. for a(pp — > H+X) relative to the SM for the total combination, 
and separately for the ee, e/j,, \x\i channels and different Higgs boson masses. 



M H (GeV) 


100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 


Exp. 
Obs. 


all: 
all: 


11.8 11.1 8.84 6.73 4.70 3.36 2.64 2.15 1.88 1.56 1.32 1.13 0.82 0.76 0.94 1.10 1.34 1.69 2.11 2.52 2.91 
17.2 19.1 13.9 8.85 5.58 4.10 2.88 2.99 2.50 2.17 1.73 1.24 0.96 0.74 0.84 0.89 1.20 1.40 2.20 2.71 2.48 


Exp. 
Obs. 


ee. 

ee 


13.6 14.2 13.5 13.0 10.1 7.21 5.56 4.32 3.86 3.22 2.77 2.30 1.76 1.64 1.94 2.24 2.72 3.44 4.12 4.90 5.65 
20.6 23.6 16.9 15.1 8.41 6.25 4.67 5.05 3.80 3.98 3.13 2.70 2.11 1.79 2.02 2.31 2.68 3.21 5.53 5.76 5.76 


Exp. 
Obs. 


e\i 


42.4 27.9 15.9 10.0 6.55 4.65 3.63 2.97 2.48 2.02 1.77 1.50 1.11 1.06 1.28 1.51 1.75 2.27 2.80 3.30 3.75 
28.1 20.0 12.5 7.95 6.23 4.75 3.39 2.88 2.61 2.14 1.75 1.28 1.00 0.86 1.17 1.33 1.79 2.16 2.68 3.31 3.20 



Exp. fifx 39.4 29.0 20.6 13.9 8.76 6.25 4.72 3.95 3.39 3.09 2.59 2.31 1.74 1.61 1.95 2.30 2.89 3.49 4.40 5.05 5.83 
Obs. nn 66.8 68.2 53.1 27.4 16.8 10.1 8.22 7.65 6.77 5.53 5.01 3.95 3.01 2.31 2.53 2.79 3.51 4.14 5.85 7.45 7.92 



TABLE V: Expected and observed upper limits at 95% CL for a(gg — > H ) x BR(H —¥ WW) in pb, for the combination of the 
ee, e/i, pp channels and different Higgs boson masses. 



M H (GeV) 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


175 


Exp. all: 
Obs. all: 


1.97 
2.55 


1.37 
2.53 


1.11 
1.60 


0.97 
1.21 


0.91 
1.15 


0.85 
0.98 


0.82 
1.00 


0.80 
1.21 


0.75 
0.98 


0.70 
0.88 


0.61 
0.69 


0.51 
0.59 


0.38 
0.42 


0.32 
0.33 


0.37 
0.34 


0.40 
0.36 


M H (GeV) 


180 


185 


190 


195 


200 


210 


220 


230 


240 


250 


260 


270 


280 


290 


300 




Exp. all: 
Obs. all: 


0.41 
0.34 


0.45 
0.34 


0.48 
0.48 


0.50 
0.44 


0.52 
0.47 


0.56 
0.61 


0.52 
0.64 


0.51 
0.65 


0.47 
0.71 


0.44 
0.67 


0.43 
0.63 


0.43 
0.72 


0.42 
0.67 


0.42 
0.70 


0.40 
0.86 





TABLE VI: Assuming fermiophobic couplings, expected and observed upper limits for cr(pp — > H + X) relative to the fermio- 
phobic Higgs expected yields for the combination of the ee, e/i, pp channels and different Higgs boson masses. 

Mh (GeV) 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 
Exp. all: 1.53 1.85 1.85 1.91 2.03 2.15 2.29 2.37 2.53 2.71 2.81 2.75 2.70 2.62 3.01 3.41 3.78 4.58 5.40 6.10 6.60 
Obs. all: 1.97 2.33 2.85 2.59 2.23 3.14 2.96 2.42 3.05 3.15 2.61 3.24 3.16 2.70 3.23 3.12 3.89 4.76 5.23 6.90 8.70 
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FIG. 19: The WW discriminant distributions in the ee channel for (a) no jet and (b) one jet. 
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FIG. 20: The WW discriminant distributions in the efi channel for (a) no jet and (b) one jet. 
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FIG. 21: The WW discriminant distributions in the fifi channel for (a) no jet and (b) one jet. 
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FIG. 22: Distribution of the combination of WW discriminants from ee, efi, and fifi channels employed to measure the WW 
cross section, sorted as a function of signal over background ratios. The hatched bands show the total systematic uncertainty 
on the background prediction. 
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FIG. 23: Results of the pp — > WW cross section measurements in each dilepton final state and their combination. 



